Orbital magnetization and its effects in spin-chiral ferromagnetic kagome lattice in the 

general spin-coupling region 
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The orbital magnetization and its effects on the two-dimensional kagome lattice with spin 
anisotropies included in the general Hund's coupling region have been theoretically studied. The 
results show that the strength of the Hund's coupling, as well as the spin chirality, contributes to 
the orbital magnetization A4. Upon varying both these parameters, it is found that the two parts 
of A4, i.e., the conventional part M c and the Berry-phase correction part Mo, oppose each other. 
The anomalous Nernst conductivity is also calculated and a peak-valley structure as a function of 
the electron Fermi energy is obtained. 

PACS numbers: 75.30.-m, 73.43.-f, 72.15.Jf 
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I. INTRODUCTION 

Recently the geometrically frustrated electron systems 
have provided hot topics in the field of condensed matter 
physics Ferromagnetic pyrochlore i?2Mo20 7 (i?=Nd, 
Sm, Gd) is one key type of the geometrically frustrated 
systems Q , which consists of corner-sharing tetrahedrons 
and the antiferromagnetic interactions between nearest- 
neighbor spins are frustrated. It was recently pointed out 
that even the ferromagnetic interaction is frustrated, if 
the spin easy axis points to the center of the tetrahedron 
In this case, the spin chirality which originated 
from the noncoplanar spin configuration, is expected to 
affect the quantum of the electrons, especially the trans- 
verse conductivity. This mechanism is also called "Berry 
phase contribution" because a non- vanishing spin chiral- 
ity is associated with a non-vanishing spin Berry phase 
for conduction electrons. To interpret the transport ex- 
periments on ferromagnetic pyrochlore 0, @|, Ohgushi 
et al. Q studied the Hall effect in a two-dimensional 
(2D) kagome lattice, which is the cross section of the 
porochlore lattice perpendicular to the (1,1,1) direction 
They obtained that if the chiral spin state is real- 
ized, the system can show a quantized Hall effect. In 
their model an important limit is used, which is that 
the electron conduction spins arc colinear with lattice 
(ions) spins. This limit is called that "the strong (or in- 
finite) Hund's coupling limit" . However, detailed experi- 
ments on the pyrochlores Q show that the chiral mecha- 
nism alone can not explain the anomalous transport phe- 
nomena in these systems. To explain these experiments, 
Taillefumier et al. [s[ studied the same lattice in a gen- 
eral case, which extrapolates the strong and weak Hund's 
coupling regions. They found that the spin Berry phase 
contribution does not depend only on the spin chirality 
but also on the strength of the local Hund's coupling. 

On the other hand, the orbital magnetism of Bloch 
electrons has been attracted renewed interest, due to the 
recent recognition [l^, [HI [HJ that the Berry phase ef- 
fect plays an important role on orbital magnetism as well 
as on the Hall conductivity. The Berry phase effect on 
orbital magnetism was until now partially presented by 



very few studies [H, [H, [3, [H, EH]- Due to its basic 
importance in understanding the magnetism and trans- 
port features of the materials, obviously, more work are 
needed in exploiting the Berry phase effect on the prop- 
erties of the orbital magnetization (OM) in various kinds 
of realistic physical systems. 

In this paper we extend the study of the OM to the 
ferromagnetic pyrochlore systems; more specially, we fo- 
cus our attention to the 2D kagome lattice with spin 
anisotropies and Hund's coupling included. It is found 
that the two parts in OM (see Sec. II), i.e., the conven- 
tional part M c and the Berry-phase correction part Mq, 
oppose each other. In particular, the OM displays fully 
different behaviors in metallic and insulating regions due 
to the different roles M c and Mq play in these two re- 
gions. Moreover, similar to the role the Hund's coupling 
plays in determining the anomalous Hall conductivity as 
observed in the above-mentioned experimental Q and 
theoretical 0] works, we find that the OM is also impor- 
tantly affected by the Hund's coupling. In particular, in 
the weak coupling case that the Mott gap bewteen the 
upper and lower Hurbard bands is overcome by the elec- 
tron kinetic energy, we show that the OM exhibits com- 
plex behaviors when scanning the Fermi energy through 
the whole series of occupied bands. Furthermore, by us- 
ing the obtained values of the OM we also calculate the 
anomalous Nernst conductivity, which is featured by a 
complicated peak- valley pattern as a function of the elec- 
tron Fermi energy. 



II. PRELIMINARIES 

Before studying the OM of the 2D kagome lattice, we 
simply review the general multiband formula for finite- 
temperature OM in the semiclassical picture of Bloch 
electrons. In the semiclassical picture [l7l. [Lsj. the Bloch 
electron for the nth band is treated as a wave packet 
|tc n (r c , k c )) with its center (r c ,k c ) in the phase space. 
The orbital magnetic moment characterizes the rotation 
of the wave packet around its centroid and is given by 
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where (— e) is the 
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charge of the electron and v is the velocity operator. By 
writing the wave packet in terms of the Bloch state, one 
obtains (k c is abbreviated as k) 

m„(k) = -i(e/2fi)(V k Unk| x [#k - 4k]|V k u„ k ), (1) 

where |w„k) is the periodic part of the Bloch state with 
band energy and is the crystal Hamiltonian 

acting on |w„k)- However, it was further found [T3| 
that the presence of a weak magnetic field B will result 
in a modification of the density of states in the semi- 
classical phase space, c£ 3 k — > c£ 3 k(l + eB-il n /H), where 
0„(k) = i(Vkit„k| x |VkM n k) is the Berry curvature in k- 
space. Due to this weak-field modification and the addi- 
tional thermodynamic average over Bloch bands included 
at finite temperature, the total free energy for an equi- 
librium ensemble of electrons in the weak field may be 
written as [T3| 
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where /i is the electron chemical potential, (3 = 1/ksT 



and e n k=e| l °k — m n(k) ■ B is the electron band 
in the presence of the external magnetic field 
equilibrium OM density is given by the field 
tive at fixed temperature and chemical potential 
- {dF/dB)^ T , with the result 



energy 
, The 
deriva- 
M = 



M=J2 Jd 3 km n (k)f n 



+ - J2 I d 3 k^n n (k) In 1 + e ^-^ 



13^ J h 

n 

M c + M n , 



(3) 



where f n is the local equilibrium Fermi function for nth 
band. In addition to the conventional term M c in terms 
of the orbital magnetic moment m„(k), the extra term 
Mfi in Eq. Q is a Berry phase effect and exposes a new 
topological ingredient to the orbital magnetism. Inter- 
estingly, it is this Berry phase correction that eventually 
enters the thermal transport current [T^]. At zero tem- 
perature and magnetic field the general expression ([3]) is 
reduced to [1(3] 

M=J2 d3k ( m «(k) + T^n(k) [mo - e„k]J , (4) 

where the upper limit means that the integral is over 
states with energies below the zero-temperature chemical 



potential (Fermi energy) /2q. 



III. THEORETICAL MODEL AND CHERN 
NUMBER 

Following previous works @, H, [H[, we consider the 
double-exchange ferromagnet kagome lattice schemati- 
cally shown in Fig. 1(a). The triangle is the one face 
of the tetrahedron. Here we consider a pure spin model 
with anisotropic Dzyaloshinskii-Moriya interactions on a 
kagome lattice. It consists of an umbrella of three spins 
per unit cell of the kagome lattice. Each umbrella can be 
described by the spherical coordinates of the three spins 
(tt/6,0), (57r/6,0), and (-tt/2,0), as shown in Fig. 1(b). 
The angle 9 ranges from to it. 

The tight-binding model of this 2D kagome lattice can 
be written as the following [j| 
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where tij is the hopping integral between two neighboring 
sites i and j; c\ a and Ci a are the creation and annihilation 
operators of an electron with spin a on the site i. Jo is the 
effective coupling constant to each local moment Si, and 
these moments are treated below as classical variables, 
rii is a unit vector collinear with the local moment S^. 
a are the Pauli matrices. In the following we change 
notation i — > (Ims), where (Im) labels the kagome unit 
cell and s is the site index in one unit cell. Note that 
in the infinite Hund's coupling limit, i.e., Jn — * oo, this 
system has already been discussed in Refs. 0, EUIlil- In 
this limit the two <r=T, J, bands are infinitely split and 
the model describes a fully polarized electron subject to 
a modulation of a fictitious magnetic flux. 

To diagonalize the Hamiltonian ([5]) , we need to rewrite 
it in the reciprocal space. We use the momentum repre- 
sentation of the electron operator 
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and the one-particle state |*(k))=XI so - * so -(k)7] o .(k)|0). 
Inserting |^(k)) into the Schrodinger equation 
H \^)=E\ 1 i), we can easily obtain the Hamiltonian 
in the reciprocal space H (k), which is given by 
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FIG. 1: (Color online) (a) Two dimensional spin-chiral fer- 
romagnetic kagome lattice. The dashed line represents the 
Wigner-Seitz unit cell, which contains three independent sites 
(A, B, C). (b) The umbrella structure on the triangular cell 
of the 2D kagome lattice. 





FIG. 2: (a) Energy spectrum of the 2D kagome lattice with 
the parameters as 8=tv/3, Jo=J c =4/\/7. (b) The critical 
value of Hund's coupling J c as a function of the chiral pa- 
rameter 6. 



where £=|ijj| as the energy unit and p k =2i cos (k • aj). 

ai = (-i,-^), a 2 =(l,0), and a 3 =(-±,^p) represent 
the displacements in a unit cell from A to B site, from 
B to C site, and from C to A site respectively. In this 
notation, the Brillouin zone (BZ) is a hexagon with the 
corners of k = ±(27r/3)ai, ±(27r/3)a 2 , ±(27r/3)a 3 , two 
of which are independent. 

Now let us consider the Chern number [2(3] and Hall 
conductivity of this system, which have been reported in 
Refs. [9l In the strong Hund's coupling limit, there 



is an energy gap between the two nearest-neighbor bands 
in the cases of 9 7^ 0, it. Then the Hall conductivity is a 
sum over occupied Bloch bands, 
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where the nth band Chern number is defined by 
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where A„ (k)=z(w n k| Vk^nk) is the Berry-phase connec- 
tion (vector potential) for the nth band. At finite tem- 
perature, considering the electron density distribution, 
the Hall conductivity is written as 

<J XV = ~ I ^f n z ■ [V k x A n (k)] . (10) 



BZ 



However in the general spin coupling cases, the gap be- 
tween two nearest-neighbor bands may disappear. When 
the Fermi energy lies in these two bands, becasue the 
gap vanishes, the Hall conductivity can not be written in 
the form of Eq. ([5]). In despite of this, the concept of 
the nth band Chern number and Eq. ([5]) are also useful 
when the gap between two nearest-neighbor bands does 
not disappear. 

For the general cases the energy spectrum can only 
be computed numerically, except for general 9 at high- 
symmetry points. For the finite values of Jo, as pointed 
in 0] , the splitting of the spectrum depends on two mech- 
anisms. One is that the coupling Jo separates each group 
of three bands, the other is that when switching on Jo, 
the pointlike degeneracies are lifted within each group 
of three bands. According to the properties of the M 
point of the Brillouin zone, a critical value of the Hund's 
coupling as a function of the chiral parameter 9 can be 
analytically obtained, which is given by Q , 



Jc(9) = ± 



2t 



VI + 3 cos 2 6> 



(11) 



Using Eq. (fTTj) one can distinguish between two different 
regimes depending on the value of Jo as compared to 
J c {9). In the regime where Jo > J c , the Chern numbers 
associated with each band are given by —1, 0, 1, 1, 0, 
— 1 from the lowest to the topmost band. Whereas in 
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FIG. 3: (a) The energy spectrum of the 2D kagome lattice, 
(b) The Hall conductivity a xy as a function of the chemical 
potential fi. In both figures, the chiral parameter is 8=n/3 
and the strength of the Hund's coupling Jo=2. The shaded 
areas are the energy gaps, which labeled as gap-I, -II, -III, 
and -IV from the lowest to the topmost gap. 



the regime where Jo < J c , the Chern numbers associated 
with each band are given by —1, 3, —2, —2, 3, —1 in the 
same order. We draw in Fig. 2(a) the energy spectrum 
with the spin chiral parameter #=7r/3 and Jo taking the 
critical value J c =4/ y7- Fig. 2(b) shows the critical value 
of the Hund's coupling J c as a function of the chirality 8. 



IV. THE ORBITAL MAGNETIZATION 

Now we turn to study the OM of the 2D kagome lat- 
tice in the general Hund's coupling cases. Similar to the 
Hall conductivity, the OM displays different behaviors in 
two regions which we will exhibit in turn. As examples, 
we consider two cases. The case I, in which we set the 
parameters as 9=tt/3 and Jo =2, is one typical case in 
the regime where J > J c (=4/ % /7 « 1.51). Wher cas the 
case II, in which the parameters are #=7r/3 and Jo=l, is 
another typical case in the regime where Jo < J c . These 
two cases are shown in Fig. 2(b) with solid dots. 

First we consider the case I. To more clearly investigate 
the OM of the 2D kagome lattice, we need to know the 
energy band structure of the system. So, we draw the en- 
ergy spectrum in Fig. 3(a), from which one can find that 
there are four gaps. From the lowest to the topmost gap, 
we denote these gaps as gap-I, -II, -III, and -IV. Clearly, 
in this case only the gap between bands 5 and 6 vanishes. 
Fig. 3(b) plots the Hall conductivity as a function of the 
electron Fermi energy (chemical potential). 

Figure 4(a) shows the OM (A4) as a function of the 
electron chemical potential fi. One can see that initially 
the OM rapidly decreases as the filling of the lowest band 
increases, arriving at a minimum at /i=— 3.54, a value 
corresponding to the top of the lowest band. Then, as 
the chemical potential continues to vary in the gap-I, 
the OM goes up and increases as a linear function of /i. 
This linear relationship in the insulating region can be 
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FIG. 4: (Color online) (a) The OM and (b) its two compo- 
nents M c (dashed line) and Mn (dotted line) as a function of 
the chemical potential /j,. The parameters are same as those 
in Fig. 3. 

understood by Eq. ([3]), from which one obtains 

A \A occu „ 

r n 

occu 

n 

Thus when the chemical potential varies in the gap-I, only 
the lowest band is occupied and dM.jd[i = — (e/h)C\. 
In this case, C\ = —1. Thus dM/dfi = e/h, i.e., the 
OM linearly increases with the chemical potential in the 
insulating region I, as shown in Fig. 4(a). Similarly, 
when the chemical potential increases in the gap-II, the 
OM increases linearly with fi. Since the Chern number 
of band 2 is zero, thus from Eq. (fT2)) and Fig. 4(a) one 
can see that the slope of the OM curve in the gap-II is 
same as that in the gap-I. When the chemical potential 
increases in the gap-Ill, the OM becomes zero and does 
not varies with the chemical potential. The reason is 
that the sum over the Chern numbers of the occupied 
lowest three bands is zero. From Eq. (fl"2"|) , one can see 
that the slope of the OM is independent of /i. In fact 
the gap-III is the usually called the Mott gap. With 
the chemical potential increases, when it lies in the gap- 
IV, one can find that the OM decreases linearly with (x. 
The reason is that the sum of the Chern numbers of the 
lowest four bands is 1. From Eq. (fl"2"|) . one can find 
dAi/d/i — —e/h, i.e., the OM linearly decreases with the 
chemical potential in the gap-IV as shown in Fig. 3(b). 

For further study, we show in Fig. 4(b) M c and Mq 
as a function of the chemical potential, their sum gives 
M. in Fig. 4(a). One can see that overall M c and Mq 
have opposite contributions to M, which implies that 
these two parts carry opposite-circulating currents. In 
each insulating area the conventional term M c keeps a 
constant, which is due to the fact that the upper limit of 
the fc-integral of m„(k) is invariant as the chemical po- 
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FIG. 5: (a) The energy spectrum of the 2D kagome lattice, 
(b) The Hall conductivity a xy as a function of the chemical 
potential fi. In both figures, the chiral parameter is 8=n/3 
and the strength of the Hund's coupling Jn=l. The shaded 
areas are the energy gaps, which labeled as the lower band 
and the higher band, respectively. 



tential varies in the gap. In the metallic region, however, 
since the occupied states varies with the chemical poten- 
tial, thus M c also varies with /i, resulting in a decreasing 
slope in the lowest three metallic regions and a increas- 
ing slope in the highest two metallic regions, as shown in 
Fig. 4(b). The Berry phase term Mq also displays differ- 
ent behavior between insulating and metallic regions. In 
the lowest two insulating region, Mq linearly increases 
with /x, and in the gap-IV, Mq linearly decreases with 
/i, as is expected from Eq. ©. In the metallic region, 
however, this term sensitively depends on the topologi- 
cal property of the band in which the chemical potential 
is located. For the bands with nonzero Chern number, 
one can see from Fig. 4(b) that Mq remains invariant, 
while for the band 2 (its Chern number is zero), it in- 
creases with the chemical potential fi. On the whole the 
comparison between Figs. 4(a) and 4(b) shows that the 
metallic behavior of M. is dominated by its conventional 
term M c , while in the insulating regime Mq plays a main 
role in determining the behavior of M . The behavior of 
the OM in this case is similar to that in the strong spin 
coupling limit. 

Then we study the case II. Similar to the case I, we 
firstly draw the energy spectrum in Fig. 5(a), from which 
one can find that there are only two gaps, which can be 
called the lower and the higher gap, respectively. In this 
case the gaps between bands 2 and 3, between 3 and 4 
(the Mott gap), and between 5 and 6 vanish. We also 
draw in Fig. 5(b) the corresponding Hall conductivity as 
a function of the electron Fermi energy. 

Figure 6(a) plots the OM (M.) as a function of the 
electron chemical potential fi. One can see that initially 
the OM rapidly decreases as the filling of the lower band 
increases, arriving at a minimum at [i=— 2.68, a value 
corresponding to the top of the lower band. Then, as the 
chemical potential continues to vary in the lower gap, 
the OM goes up and increases as a linear function of /i. 
This linear relationship in the insulating region can also 
be understood by Eq. (fT2"]) . Because the gap between 2 
and 3 and the Mott gap vanish, the OM displays a more 
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FIG. 6: (Color online) (a) The OM and (b) its two compo- 
nents M c (dashed line) and Mn (dotted line) as a function of 
the chemical potential /j,. The parameters are same as those 
in Fig. 5. 

complex behavior in this metallic region. The value of 
the OM oscillate versus the chemical potential fi. When 
\x reaches the bottom of the higher gap, the OM begins 
to go up again and linearly increases. Different from in 
the lower gap, the coefficient of the increasing is twice of 
that in the lower gap, which can be understood by Eq. 
(gH). In the higher gap dM/dfi=-(e/h) Y^=i C n =2e/h, 
while in the lower gap, dM/d[i=—(e/h)Ci=e/h. 

To see the different roles M c and Mq play in the metal- 
lic and insulating regions, we draw in Fig. 6(b) M c and 
Mq as a function of the chemical potential, their sum 
gives M. in Fig. 6(a). In each insulating area the con- 
ventional term M c keeps a constant and Mq linearly in- 
creases with fi in the two gaps with different linear co- 
efficients, as is expected from Eq. ([3]). The behaviors 
of OM in this case is novel, comparing with those in the 
infinite Hund's coupling case. There are two features in 
the regime Jo < J c . One is that the magnitudes of both 
two parts M c and Mq become much smaller than those 
in the regime Jo > J c - The other is that in the metal- 
lic region, both two parts rapidly changes. The reason 
is that the gaps between bands 2, 3 and 4 vanish and 
these bands form one energy band. However from Fig. 
6(b) one can see that the two parts of M. have opposite 
contributions to M, which implies that these two parts 
carry opposite-circulating currents. 

V. ANOMALOUS NERNST EFFECT 

The above discussion of M c and Mq can be transferred 
to study the ANE. The relation between the OM and 
ANE has been recently found (T^. To discuss the trans- 
port measurement, it is important to discount the contri- 
bution from the magnetization current, a point which has 
attracted much discussion in the past. Motivated by the 
argument [22l | that the magnetization current cannot be 
measured by conventional transport experiments, Xiao 
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FIG. 7: The ANE of the 2D kagome lattice at /fc B T=0.005. 
The parameters are (a) 9=n/3 and Jo=2; (b) 8=ti/3 and 
Jq — 1. The shaded areas are the energy gaps. 



et al. [12( have built up a remarkable picture that the 
conventional orbital magnetic moment M c does not con- 
tribute to the transport current, while the Berry phase 
term in Eq. @ directly enters and therefore modifies the 
intrinsic transport Hall current equation as follows 



Figure 7 shows a xy of the 2D kagome lattice as a func- 
tion of the chemical potential for &bT=0.005. One can 
see that the ANE disappears in the insulating regions, 
and when scanning the chemical potential through the 
bands, there will appear peaks and valleys. Remark- 
ably, a similar peak-valley structure was also found by 
the recent first-principles calculations in CuCr2Se4_ a; Br a ; 
compound [12j | . The ANE of this compound was recently 
measured by Lee et al. [H[ as a function of Br doping x 
which is used to change the chemical potential /i. Due to 
the scarce data available, until now the peak- valle y st ruc- 
ture of a xy revealed in Fig. 7 and in Ref. [H,[H,liB| has 
not been found in experiment, and more direct exper- 
imental results are needed for quantitative comparison 
with the theoretical results. Interestingly, the expression 
for a xy can be simplified at low temperature as the Mott 
relation 1131. 



t 2 B T da xy (fi a ) 



(15) 



Thus one can see that unlike the anomalous Hall effect 
[23|], ANE is given by the Fermi-surface contribution of 
the band structure and Berry curvature. Another unique 
feature of a xy is its linear dependence of temperature. 



JH = 



d 2 k 
(2tt)' 



/ n (r,k)fi n (k)-VxM fi (r), 

In the case of uniform temperature and chemical poten- 
tial, obviously, the second term is zero and the Hall ef- 
fect of 2D kagome lattice is featured by nonzero Chern 
number as discussed by Taillefumier et al. Q as well 
as in this paper. In the following, however, we turn to 
study another situation, where the current-driving force 
is not provided by the electric field (E=0). Instead, it 
is provided by a statistical force, i.e., the gradient of 
temperature T. In this case, Eqs. (fT5)) and © give 
the expression of intrinsic thermoelectric Hall current as 
jx = a xy(—VyT), where the anomalous Nernst conduc- 
tivity a xy is given by [l2l | 



1 e 



d 2 k 



Th^J (2^) 2 ' 
x (e„ k - fi) /„ + fcsTln ( 1 



(14) 



,-0(e»k-M) 



VI. SUMMARY 

We have theoretically studied the OM and ANE of the 
2D kagome lattice with spin anisotropies included in a 
general Hund's coupling region, as a supplement of the 
previous work [l5j in the infinite Hund's coupling limit. 
The results show that both of the strength of the Hund's 
coupling and the chirality contribute to the orbital mag- 
netization AL Upon varying both these parameters, it 
is found that the two parts of A4, i.e., the conventional 
part M c and the Berry-phase correction part Mn, oppose 
each other. We also calculate the anomalous Nernst con- 
ductivity and obtain a peak- valley structure as a function 
of the electron Fermi energy. 
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